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Simulation of repeaters for an optical-fiber guide digital-communication 
system requires the calculation of the statistical error rate of the signal. 
The calculations of these error rates are difficult because of the non- 
gaussian nature of the noise in the optical detector. In this paper, statistical 
techniques are described which are useful for simulating tails of distribu- 
tions. In particular, the importance-sampling procedure is used to modify 
the probability densities of the input values in a way that makes simula- 
tion possible. Application of this procedure gives more accurate results 
in reasonable computer times. The method is applied to the calculation 
of the error rate of a fiberguide repeater. Realistic examples are simulated. 
Results compare favorably with experimental measurements. The number 
of samples needed for simulation was reduced by five to six orders of mag- 
nitude. An alternative numerical solution to the problem is also developed. 

I. INTRODUCTION 

Traditional computation of error rates of digital repeaters is often 
done by assuming that the noise is gaussian in nature. For a linear 
equalizer, the error rate at a given sampling time can then be easily 
computed if the signal and standard deviation of the noise are given. 1 

The fiberguide repeater has essentially three kinds of noise sources. 
Two noise sources are located in the detector, while the third is a 
thermal gaussian noise generated by the input stage of the preamplifier 
(see Fig. 1). 

The detector, which is an avalanche photodiode is governed by two 
intrinsic noise processes : the primary photo current, which is a Poisson 
process, and the avalanche multiplication process which produces the 
gain in the diode. Neither of these processes is gaussian. 

The probability density functions (pdf) for the output current of 
the diode were derived by Personick 23 and Mclntyre. 4 The expressions 
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Fig. 1 — Optical repeater. 

describing these distributions are quite complex, and excessively time 
consuming for computer simulation of the repeater. An approximation 
to these probability densities was obtained by Webb, Mclntyre, and 
Conradi (wmc). 5 The approximate pdf (referred to below as the 
wmc pdf) is, according to the authors, fairly accurate in the practical 
operating range, and is easily implemented on the computer. 
The wmc pdf for the current at the output of the diode is : 






p(x) = 



1 



- (x - MY 



M i+ ^r>( i+ 



x-M 



(1) 



where 



M = (neq/T)G = mean output current 
q = charge of an electron 
T = time interval 
G = average avalanche gain 

n e = (vT/h£l)Popt = number of primary electrons 
t] = optical conversion efficiency 
Popt = average optical power in time T 
hSl = energy of a photon 
0-2 m n jG*F,(q/T)* = variance of the diode current 

X = Vn^VCF. ~ 1) 
F. = JfcefiG + [2 - (1/G)](1 - ft«n) = excess noise factor 
k e it = effective ionization ratio (assumed constant over range of 

G considered) 
& e ff and i\ are photodiode parameters. 

Equation (1) is, more precisely, the pdf of the mean diode current 
over the time interval T for a given average optical input power, 
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Fig. 2 — Probability density function for ONE and ZERO signals using wmc and 
gaussian approximation. 



P op t, and a given average avalanche gain G. It corresponds to the sum 
over all primary electrons of the product of the primary electron pdf 
with the avalanche gain pdf. This equation, therefore, combines the 
effects of both detector noise sources. 

In Fig. 2 the solid lines represent current pdfs for two signals: a 
ONE where the optical power is P opt = — 51 dBm and a ZERO 
where P opt = — 61 dBm. The other parameters are G = 80, tj = 0.667, 
/c e ff = 0.04, T = 22 ns. The broken lines represent error pdfs of the 
same signals and same diode parameters where a normal approxima- 
tion to the distribution is assumed. The asymmetry of the wmc dis- 
tributions is clearly seen. 

The noise generated by the first stage of the amplifier has to be 
added to this detector-generated noise. 

II. POSSIBLE SOLUTION TO THE PROBLEM 

To compute the error rate of the repeater, one has to know the 
signal probability distribution at the decision-making point, i.e., at 
the output of the equalizer in Fig. 1. 
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The output of the equalizer with impulse response h(t) is 



J — 



x(t)A(< - r)dr, 



where 



(2) 



x(<) = x.(t) + x T h(0 

x,(0 = signal output of the avalanche diode 
3th (0 - gaussian thermal noise generated by the amplifier. 

The problem is to find the distribution of y(t) for given distributions 
of x,(t) and xth(0 and impulse response h(t). 

Figure 3 shows the shapes of two distributions: P\(y) is the cumula- 
tive distribution of y for a ONE signal transmitted, P (y) is the prob- 
ability that the output current exceeds y for a ZERO signal trans- 
mitted. 2/o is the decision threshold of the comparator (Fig. 1). These 
distributions are 



Pi(y) - r p(y/o^E)dy 

J —CO 

PM = r p(y/ZERO)dy, 

Jv 



(3) 



where p(y/ONE) and p(y /ZERO) are the respective conditional 
pdfs. It should be noted that, for practical laser sources, the ZERO 
input signal could be on the order of 10 percent of the ONE signal, 
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Fig. 3— Distributions for ZERO and ONE signals. 
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since the laser is not turned off completely. In addition, there is the 
dark current of the detector diode present. The error rate of a digital 
repeater is defined as 

P e = aP 1 (y ) + fePo(yo). (4) 

(The factors a and b are determined by the a priori distribution of 
the data. For unbiased data, a = b = 0.5.) 

One of the features that has to be determined is the optimum place- 
ment of the decision threshold y 0P u that is, where P e is minimum. 

A general analytical solution to this problem is not possible ; there- 
fore, the problem has to be approximated numerically or simulated 
with the aid of a computer. Analysis of the system using a gaussian 
approximation to the detector noise sources was considered by S. D. 
Personick. 6 This analysis, however, predicts optimum decision thresh- 
olds considerably different from actual measured values. 7 Another 
solution to this problem is to replace the exact solution by a simpler 
solvable formulation that gives an upper-bound solution. This was 
done by Personick 3,8 using the Chernoff bounds. Finding upper 
bounds is a very effective technique since it gives a degree of confidence. 
Nevertheless, a method that gives more accurate results is needed 
to estimate how close the bounds are to the real solutions for this class 
of problems. 

The problem can be solved either by the numerical solution of eq. 
(3) or by statistical Monte Carlo simulation using importance sam- 
pling. Both solutions were attempted and are described below. Particu- 
lar emphasis was placed on the statistical simulation, since numerical 
solutions require excessive amounts of computer time. 

III. DISCRETE STATISTICAL ANALYSIS (DIRECT COMPUTATION OF 
OUTPUT PDF) 

The output of equalizer 7/(0 is given in eq. (2). Since we are interested 
in the output at the sampling time t , it is not necessary to compute 
the convolution for every t. The value of / is selected to maximize 

y(t ): 

y(t ) = f V x( T )h(t - t)cLt (5) 

J— 00 

since 

x{t) - Xt (t) +xth(«) 

y(to) = / x,(T)h(t - t)<1t + / z T h(t)M<o - t)<1t 

J —00 J —00 

= yM + ytaito). (6) 

The problem is to find the probability distribution of y(t ). 
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Fig. 4 — Input signal and impulse response of the equalizer. 

The probability distribution of z/th(<o) is gaussian since s T h is a 
thermally generated noise. The standard deviation of ?/th is 



ffTH 



-4 



4kTB t 

Rth 



(7) 



where B w is the effective noise bandwidth of h(t) and /2th is the 
thermal resistance. 

The effective noise bandwidth is denned as 



B «-T^f |H(/)W ' 



(8) 



where H(f) is the one-sided Fourier transform of h(t) and \H Q \ is 
the maximum absolute value of H(f). If the impulse response is 
approximated by a raised cosine of duration t, B w = 3/4t. 

Typical waveforms for x.{t) and h(t) are shown in Fig. 4. The dis- 
crete approximation of the signal y.(to) is: 



y. = ff E Xnhx-n, 



(9) 



where T is the duration of the signal and iV is the number of samples 
used in the approximation. The samples x n are independent and have 
a distribution of the form of eq. (1) for a time interval AT = T/N. 

3.1 Numerical computation ot p(y) 

If the random variables x„ are independent, then the density of 
their sum equals the convolution of their respective densities 

••(*)]■ (10) 

750 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1976 



Since the thermal noise is also additive and independent, the total 
density of the signal is 

p(y) = p(y.)*p(yTn). (li) 

3.2 Fourier transform approach 

The simplest and fastest way to compute p(y) is to transform all 
partial densities using the fft, multiply them all, and take an inverse 
fft: 

'»' - nswTO " ( ? HKafe)] 



) 



■Flv(y?My (12) 

Unfortunately, the Fourier transform is computed as a sum of positive 
and negative numbers and, therefore, suffers from large numerical 
errors. Since we are interested in the tails of the distribution (error 
rate of 10 -9 ), we conclude that the numerical errors will be excessive. 

3.3 Direct convolution 

Direct convolution produces less truncation error, because densities 
are positive functions. But direct convolution requires an excessive 
amount of arithmetic operations, e.g., if we have n distributions each 
having M points, the number of operations (multiplication and addi- 
tion) will be 

N = 5J*. (13) 

For M = 100 and n = 2, N = 10 4 which is a manageable number; 
but for M = 100 and n = 10, N = 3.5 X 10 7 and the number of 
operations becomes excessive. A program using direct convolution 
was written for n = 2. A density for the average power of the input 
pulses was convolved with the density of the thermal noise. This can 
be regarded as a response of an integrate-and-dump equalizer (capaci- 
tive input). 

The results of this convolution are shown in Fig. 5. The solid lines 
are the distributions for a ONE and ZERO being transmitted. The 
average power of the ONE in this case is — 51 dBm and of the ZERO 
-61 dBm, with T = 22 ns, « TH = 2.3 kilohms, and B w = 31 MHz, 
which corresponds to a ctth = 14.68 nA. 

The decision threshold 6 is usually expressed in percent of the mean 
range of the output current as : 

e = l°~ 1° 100 [percent], (14) 

I i ~ i o 
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Fig. 5— Distributions for ONE and ZERO signals with thermal noise ; wmc and 
gaussian approximations. 

where 

Fi = mean output current for ONE 
Y o = mean output current for ZERO 
Yo = threshold current. 

As seen from the graph, the optimum decision threshold for this 
case is 42.3 percent and the minimum error rate is 9.86 X 10~ u . 
These results were also observed in an experimental setup. 7 The broken 
lines in Fig. 5 represent the same input but with the assumption that 
the probability density of the detector current is normal. As seen from 
the graph, the resultant threshold is considerably different from the 
one obtained with the wmc distributions. 

IV. STATISTICAL (MONTE CARLO) ANALYSIS 

Monte Carlo (mc) analysis has proven over the years to be a very 
effective method for analysis of systems whose parameters or inputs 
have known statistical distributions. 9 
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However, straightforward application of mc analysis for calculation 
of error rates is impractical. The reason is that the error rates of in- 
terest are in the order of 10 -6 to 10~ 9 . To have some confidence in the 
results, 10 to 100 samples have to be in that region of low probability. 
This will require on the average 10 11 samples, which is impractical. 

To overcome this shortcoming, we introduced a modified sampling 
method called importance sampling. 10 Conceptually, the idea is simple. 
If the regions which contribute to the result are known, modify the 
distributions in such a way that more samples are taken from the 
important regions. 

Let us consider the situation where the problem is to estimate the 
value of 

S = r g(x)p(x)dx, (15) 

J— oo 

where p{x) is a pdf. The importance-sampling procedure consists in 
introducing another pdf p*(x), which is preferable for sampling pur- 
poses. Then 

= / g*(x)p*(x)dx. 
£ is then estimated from 

where x, is picked from the pdf p*(xi). 

4.1 Defining the Important regions for error-rate computation 

The problem is to define the regions where the errors come from. 
In Fig. 6a is shown the input waveform for a ONE, while Fig. 6b 
depicts the impulse response of the equalizer and 6c is the pdf of the 
diode current for one time interval AT. 

Since the input signal is always positive, it is obvious that for posi- 
tive h„ most of the errors, that is, outputs that are below the threshold 
0, will come from region A in the pdf of Fig. 6c. Conversely, if h n is 
negative, most of the errors will come from region C. Exactly the 
reverse is true for transmitted ZERO signals, since the signal is con- 
sidered erroneous if it exceeds the threshold 6. 

A very similar argument can be constructed for thermal noise, 
where the "important" regions are the positive tail for a ZERO trans- 
mitted and the negative tail for a ONE transmitted. 
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Fig. 6 — (a) Optical input signal, (b) Impulse response of the equalizer, (c) pdf 
of average detector current for AT. 



4.2 Biasing of the PDF 

To implement the importance-sampling technique, the distributions 
have to be modified so that the important regions are emphasized. 
The pdfs in the random-number generator that we use are stored in 
tabular form, as shown in Fig. 7a. The region in question can be 
biased by a constant multiplier as in Fig. 7b, where p* = pB for 
o ^ x ^ b and where the bias B is constant. In Fig. 7c, the pdf is 
biased such that it is quasi-uniform in the region prescribed. The bias 
Bi is constant in intervals between two adjacent points in the table. 
For each interval (x,_i, a\), this bias is 



Bi = 



<xp, 



P< 



(16) 



where p, is the value of the pdf at x, and a is a coefficient specified by 
the user. The biasing and normalization is done automatically by the 
program. The output of the random-number generator consists of 
two values: a random number RV and its associated bias BV. BV is 
the bias of the interval that RV was picked from. 
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Fig. 7 — Probability density functions, (a) Unbiased, (b) Constant bias, (c) Bias 
for quasi-uniform pdf. 



4.3 Computing the histogram 
4.3.1 The usual histogram 

The standard output of a Monte Carlo analysis is the histogram. 
Usually the histogram is constructed by subdividing the output vari- 
able into M bins, Fig. 8. 

The statistical analysis program counts the number of outputs which 
fall into each bin. For a large number of samples (N —>«>), the histo- 
gram approximates the pdf. The probability p, that the output y will 
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Fig. 8 — Histogram, 
be in the interval y,_i S y = y< is 



where 



Pi = lmi m ' 

N-> 00 J\ 



n,- = number of samples in bin i 
N = total number of samples. 



(17) 



The accuracy of estimation of pi can be established by 11 



-?-W 



(1-p) 



(18) 



p \ JVp ' 

where ±e is the 95.4-percent confidence limit of p i} p is the desired 
probability, and <r, is the standard deviation of p,-. 

4.3.2 Histogram tor Importance sampling 

Let us assume we have M independent random variables Xi, • • ■ , 

y = f(xi, •• •, x m , • • •, x M ), 

p*(x m ) is the biased pdf of the random variable x m . How do we com- 
pute the unbiased histogram for y. 

If a sample x m i is selected from p*(x m ), its probability to be selected 
is increased by the bias B mi . Therefore, the weight of the sample x m i is 

1 



W mi = 



B, 



Since y, is determined by M independent variables x m the probability 
of yi, to be selected from these regions, is increased by the product of 
all the biases of x m . Therefore, the weight of yt is Wi = 1/(11 B mi ). The 

m 

average weight of all the n, samples that fell into bin j is given by 

w "''hk v "'"h%XwB:)i (19) 
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The probability that the unbiased output will be y,-i ^ y ^ y, is 

Pi = w« v y -^ = w.v/P*. (20) 

Since this method increases the probability p in eq. (18), the error 
t is reduced considerably. 



- 2 -'^ - 2 Vt- (21) 



For error-rate computation p « 1 , 



~ = 2*fe . (22) 



ViVp* \Np' 

and since p* is increased in the interesting areas by l/to aT /, the con- 
fidence limit is decreased on the average by 

7 = — = V^~ (23) 

where 

eis = Monte Carlo analysis with importance sampling 
ec = conventional Monte Carlo analysis. 

4.4 Example 

A simple example to illustrate the principle of importance sampling 
is the problem of calculating the probability of obtaining a total of 
three when one tosses two dice. Each die has six faces labeled from 
one to six and each face has the probability p = 1/6 of being on top. 

The problem, of course, can be solved analytically. Any particular 
combination of the dice has the probability of 1/6-1/6 = 1/36 of 
occurring. Since there are two combinations which make three (one- 
two and two-one), the probability of getting a three in a random toss 
of the dice is 2/36 = 1/18. 

If the probabilities of a one and a two are biased so that they occur 
twice as often as usual, that is, p = 1/3 rather than 1/6, then the 
probability of getting a three will be four times as great as in the un- 
biased case. That is p = 2/9 instead of 1/18. The weight for this 
example is 

w, = 1/(B 1 B 3 ) = 1/4 



and the probability of obtaining a 3 is being approximated by 

1 >H 

In' 
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£ 1 /ts 
P = 7 at 1 



The error in eq. (18) is reduced by slightly more than a factor of two, 
since for 

8.24 



and for 



V = 1/18, e = 

V = 2/9, e = 



Viv 

3.74 

Vat ' 



V. ERROR-RATE ANALYSIS 

The importance-sampling method described above was applied to 
the error-rate analysis of the fiberguide repeater with an avalanche 
diode detector. 

To demonstrate the analysis procedure, we will analyze one 
specific case. Let us assume that the input to the photodiode is as 
shown in Fig. 9. The ONE is a square pulse of —51 dBm power of 
22 ns duration. ZERO is 10 percent of the ONE power or —61 dBm. 
The avalanche diode data are 

average avalanche gain = 80.0 
ionization ratio k = 0.04 

conversion efficiency 17 = 0.667. 

The equalizer impulse response h(t) is as in Fig. 10 and the thermal 
input resistance is Rrn = 2.3 kilohms. 

The pulses (ONE and ZERO) are subdivided into nine intervals 
[eq. (9)2, since this resolution is considered adequate for the impulse 
response in Fig. 10. Nine distributions* of the form in eq. (1) are com- 
puted for the ONE and ZERO pulses. 



P 0PT IN dBm 



T = 22 ns 
Fig. 9 — Optical input. 



* In this case the distributions are identical since the pulse is rectangular. 
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Fig. 10 — Impulse response. 

The mean for ONE is M x = 400.3 nA 

The standard deviation for ONE is a x = 102.66 nA 

The mean for ZERO is M = 40.03 nA 

The standard deviation for ZERO is a = 32.46 nA 

The effective noise bandwidth []eq. (8)] for the thermal noise is 
B w = 31 MHz. The standard deviation is computed from eq. (7) : 

(7th = 14.68 nA. 

From the above considerations, the output for N intervals is 



T N f N 

2/ = Tr L X„h N -„ + -rr £ /l„X T H, 

N „-i Jy „=i 



(24) 



where the x„ and xth are samples from the respective pdfs and h's 
are discrete values obtained from h (t) . The impulse response is normal- 
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ized to produce unity gain 

J 1 N 

K= ^=_ (25) 

£ /in 

n = l 

AT 

V* = £ On + ZTH- 
ji = 1 

In addition, the above case has been analyzed for N = 1. This cor- 
responds to the integrate-and-dump equalizer which was analyzed 
numerically in Section 3.3. 

5.1 Computation of distributions and biasing 

The calculations of distributions for the diode output are straight- 
forward using the wmc distributions of eq. (1). The biasing is a more 
complex procedure and sometimes requires a few trials. The idea is 
to get a reasonable amount of samples in the "important" region of 
the output. The pdf for ONE was computed for x„ = to 800 nA, 
and was biased according to eq. (16), 

B { = ot ^^ for to 400 nA with a - 1. 
Pi 

The pdf for ZERO was computed from to 400 nA and was biased 
from 40 nA to 400 nA with an a = 0.1. The gaussian distribution for 
the thermal noise was tabulated for ±7<r. The biased regions were 
-7<r to for the ONE signal and to 7a for the ZERO signal with 
an a = 0.5. The pdfs were stored in tabular form for 100 points. Since 
the ZERO signal has a pdf which requires more detail for small cur- 
rents (Fig. 2), the independent variable was spaced logarithmically. 
All others were linearly spaced. 

The distributions were computed by a discrete form of eq. (3) : 

Cuml(n) = £ Histl(i) 

(26) 

M 

Cum0(n) = £ HistO(i), 

i=n 

where Histl(z') and Hist0(t) are histograms of the outputs. 

To calculate the minimum error rate of eq. (4), Cuml(rc) and 
Cum0(n) were smoothed and approximated by asymptotic func- 
tions. The minimum of the sum of these functions is the minimum error 
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rate, and the current where this minimum occurs is the optimum 
threshold. 

5.2 Smoothing of the cumulative distributions 

It was shown by S. B. Weinstein 12 that for the class of "exponential- 
type" distributions the error-probability distribution has an asymptote 
in the tails of the form 

P.Cr)~ex P r-^Yl- (27) 



For v = 2 the distribution is gaussian. 

To obtain better accuracy, this asymptotic function was expanded 
to 

P e (x) = exp{ — exp[a(log x) 2 + b log x + c~\) 

or (28) 

log ( — log P e ) = a(log x) 2 + b log x + c ; 

i.e., the log log of P e is a second-order polynomial of log x. This equa- 
tion reduces to Weinstein's form, eq. (27), for a = 0. 

The second-order polynomials were then least square fit to both 
Cuml and CumO. 

5.3 Minimum of error and optimum threshold 

The probability of error, eq. (4), is defined as 

Pe(D = KP«d) + P«(D1 

= Hcxpf-expfXlog/) 2 + 6 1 (log 7) + cQ\ 

+exp{-exp[a (log7) 2 + 6 (log 7) + c ]}), (29) 

where 7 is the decision threshold current. 

To find the minimum of P,(I), we compute 7 op t from 

^1 = 0. (30) 

Pe(I)min = Pe(^opt)- 

Equation (30) is solved using the Newton-Rap hson method. The 
threshold 6 is usually expressed in percent as in eq. (14). 

VI. RESULTS 

Figure 11 shows a histogram for the ZERO signal for 10,000 samples. 
The abscissa is scaled logarithmically. In the tables on the right side 
are indicated the ordinate values for the bins, the number of samples 
in each bin, and the distribution for the bin. Below the histogram are the 
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Fig. 11— Histogram for ZERO signal. 

coefficients of the fitted polynomial R(I) = a + a x log / + a 2 (log I) 2 , 
where R(I) = logio [-logio-P eo CO]- 

Figure 12 shows the smoothed distributions for the above example 
for one section (N = 1) and 1000 samples. The stars represent nu- 
merical solutions for the same input parameters, as described in Section 
III. As seen from the data, if biasing is chosen rightly, a comparatively 
small sample size is adequate to simulate the extreme tails of the pdf. 

Figure 13 shows the smoothed distributions when h(t) is represented 
by nine sections (N = 9). The broken lines represent the distributions 
for N = 1. The sample size was 10,000 in both cases. As seen from 
the graph in this example, an integrate-and-dump equalizer analysis 
can be used as a first-order approximation to an equalizer design with 
an impulse response similar to the h{t) shown in Fig. 10. 
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The minimum error-rate and optimum threshold were 



N = 9 P emin = 3.03 X 10- 1 
N = 1 P emi „ = 6.6 X 10- 11 



= 43.61 percent 
= 42.68 percent. 



Figure 14 again shows the distribution of the above example (N = 9). 
The circled points are results obtained using Chernoff bound approxi- 
mations as obtained by J. H. Bobsin 8 using Personick's 3 method. As 
seen from the graph, the optimal threshold predicted by the Chernoff 
approximation is identical with the one obtained by the statistical 
analysis. The minimum error rate shown by the Chernoff bound is 
about two orders of magnitude larger from the statistically simulated 
value. However, the error rate obtained by the Chernoff approximation 
is an upper bound. Therefore, if the result falls within the acceptable 
design specifications, no further analysis is necessary. If, on the other 
hand, the predicted error rate is marginal, a more accurate analysis 
technique such as the one described in this paper is advisable. 

Figure 15 depicts curves for the minimum error rate P em in and 
optimum threshold as a function of the standard deviation of the 
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Fig. 12— Distributions for ONE and ZERO signals. 
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Fig. 13 — Distributions for ONE and ZERO signals using Monte Carlo method for 
1 and 9 sections. 
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Fig. 14 — Distributions for ONE and ZERO signals using Monte Carlo method and 
Chernoff bounds. 
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Fig. 15 — Minimum error rate and optimum threshold for 1 and 9 sections as a 
function of thermal noise. 



thermal noise o- T h. For this example, 

23.32 



0"TH — 



Vb 



!>A], 



TH 



where Rth is the input impedance of the equalizer in kilohms. 

The solid lines represent P e win and op t for N = 9. The broken lines 
depict the same case analyzed with N = 1. 

As seen from the graphs, the optimum threshold for both types of 
receivers is between 42 and 46 percent. Interestingly, the threshold has 
a minimum at a = 25 nA (Rth = 0.87 kilohm). 

The integrate-and-dump analysis has a lower error rate then the 
nine-section equalizer analysis in the normal operating range of the 
repeater. 

Throughout this analysis, the avalanche gain was assumed to be 
constant. In a more realistic design, the avalanche gain can be adjusted 
to optimize the error rate. This will result in a higher avalanche gain 
for higher thermal noise. 

In this example, the input signals and the equalizer impulse re- 
sponse were chosen so that there is no intersymbol interference. This 
was done to simplify the example. However, the analysis procedure 
is not limited to this case. 
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VII. CONCLUSIONS 

The statistical importance-sampling technique was used to calculate 
error rates of optical repeaters. The application of this technique sub- 
stantially reduced the number of samples needed for simulation of the 
repeater and thus made computer simulation practical. 

The importance-sampling technique is very effective and has general 
applicability in areas where tails of distributions have to be computed 
or particular sensitive areas of the statistical distribution have to be 
emphasized. 

This method cannot, however, be used automatically ; a good prior 
knowledge of the behavior of the system is necessary. To get an 
accurate statistical simulation of a section of an output, all the im- 
portant regions of input pdfs which may contribute to this section 
have to be emphasized. 
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